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Abstract. We consider the classical picture of three dimensional motion of charged particles in pulsar magneto- 
sphere. We adopt a perturbative method to solve the equation of motion, and find the trajectory of particles as 
they move along the rotating dipolar magnetic field lines. Our aim is to study the influence of rotation on the 
pulsar radio emission by considering the constrained motion of particles along the open dipolar magnetic field 
lines. We find that the rotation induces a significant curvature into the particle trajectories. Our model predicts 
the intensity on leading side dominates over that of trailing side. We expect that if there is any curvature induced 
radio emission from the region close to the magnetic axis then it must be due to the rotation induced curvature. 
Our model predicts the radius-to-frequency mapping (RFM) in the conal as well as core emissions. 
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1. Introduction 

The process behind the radio emission in pulsars still 
evades a concrete understanding. The different types of 
radiation mechanisms proposed for pulsars are only par- 
tially successful in explaining the phenomenon of radio 
emission. Out of them, the curvature emission emerges as 
a favorable model which can satisfactorily explain some 
aspects of pulsar radiation like polarization angle swing. 
Though it has unresolved issues like the bunching and co- 
herent of emission, it is still believed to be a natural and 
unavoidable emission process for pulsars (Gil et al. 2004). 
In this perspective an effort towards understanding the 
curvature emission process is worthwhile to pursue. 

A precise knowledge about the particle motion is re- 
quired for understanding the curvature emission in pul- 
sars. Over the past three decades almost all the works, 
which considered the curvature emission in pulsar mag- 
netosphere, assumed that the trajectories of particles are 
identical to the dipolar field lines along which they are 
constrained to move. The important effect, which is ne- 
glected in this approximation, is the influence of pulsar 
spin on particle trajectory. However, for an observer in 
the laboratory frame (inertial frame), the curvature ra- 
dius of the particle trajectory need not be same as that of 
the field lines along which they move. Although there are 
occasional notes about the rotation induced curvature into 
the particle trajectory (e.g., Hibschman & Arons 2001), a 
proper estimation is still lacking. 



We restrict our analysis to the region of magnctosphcrc 
where the radio emission is prominent. Some limits on 
the radio emission heights have been estimated by geo- 
metrical methods (e.g., Gil & Kijak 1993; Kijack & Gil 
1997), interstellar scintillation (e.g., Smirnova & Shishov 
1989), and aberration phase shift in pulse profiles (e.g., 
Blaskiewicz et al. 1991; hereafter BCW91). While dis- 
cussing the concept of radius-to-frequency mapping, sev- 
eral authors (Phillips 1992; Kijak & Gil 1997; Kramer 
et al. 1997) have set the constraints on the range of radio 
emission altitudes. Kijak and Gil (1997) set an upper limit 
of 10% tl for the normal pulsars, where tl = Pcj (27r) is 
the light cylinder radius, P is the pulsar period and c is the 
speed of light. Recently, Gangadhara and Gupta (2001) 
have analyzed the emission geometry of PSR 0329+54 and 
estimated the aberration-retardation (A/R) phase shift of 
pulse components by presupposing a nested conal struc- 
ture for the emission beam. They have found that the ra- 
dio emission altitude is not constant across the pulse win- 
dow but it varies in such way that inner cones including 
the core comes from lower altitude (< l%r\) compared 
to the outer ones (~ 5% rt)- In a subsequent work Gupta 
and Gangadhara (2003) have verified it for six more pul- 
sars. Based on these deductions wc surmise that an upper 
limit of 10% n, for the radio emission heights of normal 
pulsars and 20% tl for millisecond ones (e.g., Gangadhara 
2005a). 

We solve the equation of motion and analyze the dy- 
namics of charged particles moving along rotating mag- 
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nctic field lines. The equation of motion of a particle in 
the combined electric and magnetic field is given by 

dp\ a b 



dt 



= qE + F B , 



(D 



where p lab is the momentum of particle in the laboratory 
frame, q E and Fb are the electric and magnetic parts of 
the Lorentz force, and q is the particle charge. 

The mechanism of initial acceleration of particles in 
the regions of a few stellar radii has been analyzed by sev- 
eral authors (e.g., Ruderman & Sutherland 1975; Arons & 
Scharlemann 1979; Harding & Muslimov 1998; hereafter 
HM98). A simplified picture of initial acceleration, which 
one can gather from these models, may be summarized 
as follows. The huge electric fields, which are rotation- 
ally generated very close to the neutron star, accelerate 
the primary particles against gravity to acquire very high 
Lorentz factors of as 10 5 — 10'. These primary particles 
in turn emit the pair-producing 7— ray photons through 
curvature radiation. These photons in turn produce the 
secondary electron-positron (e + ,e~) pair plasma, which 
further develops into pair formation fronts. By developing 
a highly sophisticated model for the initial acceleration of 
charged particles HM98 have shown that secondary par- 
ticles can screen the accelerating electric fields (Eh = 0) 
above a distance of about 1-2 stellar radii from the neutron 
star surface, except in the regions of slot gaps (Muslimov 
& Harding 2004). Thus we take that the parallel electric 
field E|| = in the radio emission region. 

Ruderman and Sutherland (1975, hereafter RS75) 
have given an estimate for the Lorentz factor (7 as 800) 
of the secondary particles when they come out of the 
region of huge initial acceleration. Recently, by solving 
the viewing geometry in dipolar dipolar magnetic field 
lines, Gangadhara (2004) has deduced Lorentz factor of 
(7 as 300) secondary pair plasma. So, we take that the 
particles, which escape from this region, have Lorentz 
factor 7 as 100-1000. The charged particles, which enter 
the region of .Em =0, experience an acceleration primarily 
because of curvature in their trajectory, and hence emit 
curvature radiation. 

The multipolar components of pulsar magnetic field 
are expected to be limited to very low altitudes close to 
the stellar surface and vanish at higher altitudes leaving 
a dipolar configuration. Observations indicated that the 
magnetic field in the radio region is more likely to be dipo- 
lar (Xilouris et al. 1996; Kijack & Gil 1997). An estimate 
for the azimuthal distortion of magnetic field due to mag- 
netic field line sweep back is given by Dykes and Harding 
(2004): 



xy I x 

o—^- cos a + 3— — 1 I sma 



,(2) 



where r n = r/r-^. x = r cost// sin#', y — r sin0' sin#' 
and z = r cos 9' . The angles <fi' and 9' are the azimuth and 
colatitude of the emission spot in a spherical co-ordinate 
system centered on the rotation axis. The radial distance 



r is measured from the neutron star center to the emission 
spot. The value of 5<?Wsb is found to be much less than 
the A/R effects in the region of our interest (Gangadhara 
2005; hereafter G05). Thus we take dipolar magnetic field 
is represented by 



B 



3(f • m)f — m 



(3) 



where m is the magnetic dipole moment of the star, and 
the field lines are governed by 



r c sin 



(4) 



where 9 is the magnetic colatitude and r c is the field line 
constant. Large values of r c represent the field lines with 
larger radius of curvature and vice versa. 

The drift velocity v cg of charged particles across the 
field lines due to curvature and gradient in the magnetic 
field (Jackson 1975) is given by 



2v 



2Ruj b 



Rx B 

RB 



(5) 



where v\\ and v± arc the parallel and perpendicular com- 
ponents of particle velocity with respect to the mag- 
netic field B, R is the radius of curvature of field line, 
wb = qB/(mc) is the cyclotron frequency, m and q are 
the mass and charge of particle, respectively. Based on the 
estimate given in BCW91, we infer vh ~ c and v± <C c. 
For the typical pulsar parameters we find v cg -C W11, and 
hence we neglect the drift motion. 

Machabeli & Rogava (1994) by considering the parti- 
cles moving freely along an infinitely long and rigidly ro- 
tating straight tube, derived an expression for the particle 
trajectory. Reiger & Mannhiem (2000) have also discussed 
the particle acceleration along the rotating straight mag- 
netic field lines in AGN, by assuming the angular velocity 
of particles is same as that of AGN. Further, Rogava et al. 
(2003) have considered the free motion of particles within 
an arbitrarily curved and infinitely long rotating cylinder. 

In our previous work (Thomas & Gangadhara 2005; 
hereafter TG05), we considered the two dimensional (2D) 
geometry and analyzed the dynamics of a charged particle. 
We approximated the projected field lines on the equato- 
rial plane to be straight lines. However, in the present 
work we relax this assumption by considering the three 
dimensional (3D) motion of charged particles in the ro- 
tating dipolar field lines. In §2 we present the solution to 
the equation of motion, and its solution. We discuss the 
possible implications of our model in §3. Concurrently, we 
also find a numerical solution to the radial position of par- 
ticle as a function of time in §3.3, and conclude in §4. The 
detailed derivations are provided in the appendices ( A-E) . 

2. Equation of motion and its solution 

We derive and analyse the equation of motion of charged 
particle moving along the rotating magnetic field lines in 
the radio emission region. The derivation of equation of 
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motion and its solution are given in detail in Appcndiccs- 
B and C. Using Eq. (|C-20[) . we rewrite the equation of 
motion Eq. (|C-28[) as: 



2 n 2 m rr 2 + {dnlJdt) rr 2 



c 2 (l-Q 2 m r 2 /c 2 ) 



CC r = 



where 



(6) 



(7) 



r = dr/dt and r = d 2 r/dt 2 . Equation ([6]) is highly nonlin- 
ear as it contain the terms in f 4 , and the exact analytical 
solution seems too difficult to achieve. So, we opt for a 
perturbative method to find solution. Let e = r^/r e be 
the perturbation parameter. We are limiting the consid- 
eration of particle motion to the open field lines for which 
r c > 5 tl for the reasons explained in Appendix-C. After 
substituting Eqs. (|D-1[) - (|D-6|) into Eq. the coefficients 
of e° give the zeroth order equation: 



ro 



c 2 (l-n 2 nQ r 2 /c 2 ) 
and those of e 1 give the first order equation: 



4n 2 m0 r r t ( 2n 2 m0 f 2 

' 1 I q , „ 9 9 / q \ ' 1 



(8) 



c 2 (l-^ io7 .2 /c2) Vc^l-f^/c 2 ) 



c 2 (l-n 2 m0 r 2 /c 2 ) 

2n 2 n n 2 ifM 

mO ml 



(9) 



In our previous work (TG05), we have already found 
the solution of the zeroth order equation: 



r = — — cn(A - ft m0 t) , 
"roO 

and to the first order equation, the solution is 



i V2 k I y\ ^ 
n = -yi / dt + y 2 / at , 

w w 



(10) 



(11) 



where the expressions of the parameters yi, yi, w and k 
are given in Appendix-D (see Eq. D 21). 



3. Discussion 

Based on our solution to the radial position of the charged 
particle, we discuss some of the features of particles tra- 
jectory and their radius of curvature relevant in the radio 
emission region of pulsars. 



3.1. Particle Trajectory 

Consider a Cartesian coordinate system (x, y, z), whose 
origin is fixed at the neutron star center and z-axis taken 
along the rotation axis (see Fig. 1). Using Eqs. (|D-1|) . 
(HSU), (|T>20|) . and (jA~T|) and (|A~4)) given in Appendix- A, 
we find the trajectory of particles: 

r = ?-(i)(siii 9' cos 4> p e x + sin^' sm<f> p e y + cos 9' e z ) , (12) 

where (e x , e y , e z ) are the unit vectors along the (x, y, z) 
axes, respectively. 

In Fig. 2 we have plotted the trajectory of particles 
moving along the rotating field lines specified by field line 
constant r c = 5 ?*l and magnetic azimuth 4> = ±90° in 
the case of an orthogonal rotator. We have plotted the 
sequential locations of a rotating field line with time step 
of 0.0055 sec in an inertial frame. The upper-half of the 
diagram shows the leading side and lower half shows the 
trailing. The intersection points between the trajectory 
and field line mark the instantaneous positions of particle 
during its passage. The continuous solid line, which joins 
those points, define the trajectory of particle in the labo- 
ratory frame (inertial frame) . It is clear that they bend to- 
wards the direction of rotation (counter-clockwise). Since 
the field lines on trailing side bend in the opposite direc- 
tion to rotation, we see an asymmetry in the trajectories 
between the leading and the trailing sides. In the absence 
of rotation, trajectories are expected to symmetrically di- 
verge from the magnetic axis as the field lines do, but in 
the presence of rotation they become asymmetric. It is 
to be remarked that even for a slight addition of trans- 
verse velocity, because of rotation, to the radial velocity 
(u|| ~ c) the change that it brings in the radius of curva- 
ture is significant. An extreme example would be the case 
of a particle moving along the rotating magnetic axis. At 
an altitude of ~ 0.1 Tl the value of radius of curvature 
p = oo if no rotation, but with rotation velocity (<0.1c) 
it would be p ~ rL/(2 sina). 

In Fig. 3, the trajectories of particles moving along the 
magnetic axis (6 — 0°, & — a and 4> = 0°) are plotted for 
the different cases of inclination angle a. We infer, the tra- 
jectories are curved towards the direction of rotation, and 
the radius of curvature of the trajectory increases with the 
decreasing a, as the rotation effects are more significant 
in the case of orthogonal rotators. It demonstrates the 
influence of rotation on the trajectory of particles in lab- 
oratory frame. Even though the magnetic axis is straight, 
the trajectory of particles are not as they travel in time 
from slowly rotating region (polar cap) to fast rotating re- 
gion (light cylinder) . We note that the relativistic particles 
moving in such curved trajectories emit the curvature ra- 
diation. To further illustrate the influence of rotation on 
particle trajectories, we have plotted them on the equa- 
torial plane (see Fig. 4) of a pulsar with a = 90° and 
P = 1 sec. We have selected two open field lines r c = 5r^: 
one on the leading side (<j> = 90°) and other on the trail- 
ing side (4> = —90°). For reference we have also plotted 
the trajectory of particle moving along the magnetic axis 
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Fig. 1. The geometry showing the motion of a particle along a 
rotating magnetic field line (thick curve PO) in the laboratory 
frame. The phase difference A<f> is an angle between the pro- 
jected magnetic axis AO and the projected radial vector OB on 
the XY-plane. The angle IXOA— Qt is the phase of projected 
magnetic axis with respect to X-axis of a fixed Cartesian co- 
ordinate system-XYZ whose Z-axis is parallel to the rotation 
axis Cl. 



(4> = 0°). Based on Figs. 2, 3 and 4 we make the following 
observations in the laboratory frame: (1) the trajectories 
arc deflected towards the direction of rotation with respect 
to the field lines, and the magnitude of deflection increases 
with the radial distance from the rotation axis. (2) The 
magnitude of deflection on leading side is larger than the 
trailing side, i.e., an asymmetric deflection seen. This is 
because of the fact that the dipolar field lines bend toward 
the direction of rotation on leading side whereas on trail- 
ing side they bend in the opposite directions. Hence we 
expect the rotation induced curvature emissions on lead- 
ing side to dominate over the trailing side. 

3.2. Radius of Curvature of Particle Trajectory 

In a corotating frame particles are expected to move with 
relativistic velocities V|| ~ c along the field lines with 
Lorcntz factor 7 ~ 100 — 1000. But for an observer in the 
Laboratory frame particles move rclativistically along the 
field lines as well as rotating with the velocity v rot = fixr. 
The emitted radiation will be beamed in the direction of 
total velocity: vi a b = vm + v rot . Since rotation velocity is 
a function of radial distance r from the rotation axis, it 
increases with altitude. It is quite small at the polar cap 
and approaches the speed of light (c) near the light cylin- 
der. In the radio emission zone (r ~ 200 to 2000 Km) of 
normal pulsars, the rotation velocity is about 2 to 10 per 
cent of c for the spin periods of 1 to 0.1 sec. Further, in 
millisecond pulsars it is still larger. Therefore, with an in- 




0.05 0.1 0.15 0.2 

x/r L 

Fig. 2. A diagram showing how the rotating field lines in- 
fluence the particle trajectory on the leading (L) and trailing 
(T) sides of the magnetic axis. The several lines on the upper- 
half of the diagram represent the sequential rotated positions 
of a single field line with a time step of 0.0055 sec. Similarly, 
those on the lower-half the figure mark the sequential loca- 
tion of field lines on trailing side. The intersections (marked 
with bullets) between the field line and trajectory specify the 
instantaneous location of particle in the observer frame (iner- 
tial frame). These trajectories are computed by considering the 
field lines specified by r c = 5 tl and (f> = ±90° on the equato- 
rial plane of an orthogonal rotator having spin period of 1 sec. 




Fig. 3. The trajectories of particle moving along the magnetic 
axis at different inclination angle 01. The azimuthal distortion 
due to the magnetic field sweep back is neglected in this plot. 
Used pulser period P = 1 sec and = 0°. The curve labeled 
with a = 90° lies in the xy-plane, and all other lie above it. 



creasing r the trajectory of particles, which corotate with 
the magnetic field lines, become more and more curved to- 
wards the direction of rotation. Even if we take a straight 
field line, the trajectory of particle becomes curved for an 
observer in the laboratory frame, as the particle gets dif- 
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Fig. 4. The projected trajectories of particles moving along 
the field lines with <j> = 90° on leading side, and (f> — —90° on 
trailing side, on the equatorial plane. The solid lines represent 
the magnetic field lines while the dotted ones mark the particle 
trajectories. Chosen a = 90° and P = 1 sec. 
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-0.05 0.05 0.1 0.15 0.2 0.25 

x/r L 

Fig. 5. The curvature radii for the trajectories L and T on 
the leading and trailing sides (see Fig EJ . At x — 0.05 tl the 
computed values of curvature radii are p = 0.2 tl and 0.4 tl 
with local centers at Ol and Ot for the leading and trailing 
side trajectories, respectively. 



ferentially shifted towards the direction of rotation. This 
introduces a real physical curvature into the trajectory of 
particle in the laboratory frame. Though at any particular 
instant an inertial observer can see only a small segment 
of the trajectory, as the time progresses he will be able to 
see the full trajectory due to pulsar rotation. 

The curvature radius of particle trajectory can be es- 
timated using the following expression 



where v and a are the particle velocity and accelera- 
tion, respectively. In Fig. 5, we have plotted the trajec- 
tory of particles and their curvature radii on leading and 
trailing sides as fractions of light cylinder radius n,. At 
x = 0.05 tl, we have marked their local center of cur- 
vature and curvature radii. We find curvature radius on 
leading side (0.2 rx) becomes smaller than that on trailing 
side (0.4 tl). The "projected acceleration" can be given 
as a pro j = v 2 1 p ?s c 2 /p, where v is the instantaneous ve- 
locity of particle and p the radius of curvature, which is 
influenced by the rotation. Hence the emissions on lead- 
ing side are expected to be stronger than the trailing side, 
in agreement with the results anticipated by Blaskiewicz 
et al. (1991) and Peyman & Gangadhara (2002). 

By considering the motion of particles along the field 
lines, which lie in the meridional plane defined by (Cl,rh), 
we have plotted the curvature radius as a function of r/r^ 
in Fig. [6l It shows that when the rotation is taken into ac- 
count the curvature radius of particle trajectory becomes 
far lesser than the intrinsic curvature of the field lines. The 
value of p approaches the limiting value « nV(2 sin a) at 
higher altitudes for the field lines closer to the magnetic 



axis, as the rotation picks up and induces curvature into 
the trajectory. The curves are roughly flat at higher alti- 
tudes (r/rx, > 0.01) because of the fact that the intrinsic 
radius of curvature of field lines increases with altitude 
but the rotation induced radius of curvature decreases as 
the rotation picks up at larger altitude. Though the in- 
herent curvature radii of these field lines are much above 
tl/(2 sin a) , it is the rotation which brings radius of curva- 
ture of particle trajectory to this limiting value of ~ t\l/2. 
It is also evident from the approximate expression of p m 
(see below) . For Fig. we considered the field lines which 
lie outside the meridional plane (0 = ±90°). It shows due 
to rotation the curvature radius on the leading side be- 
comes smaller than that on the trailing side. Since the 
field lines having large (say r e > 100) are almost straight, 
the difference in curvature radii between the trailing and 
leading sides is minimal at higher altitudes (r/r^ > 0.01). 

In most of the related works the curvature of particle 
trajectory has been assumed to be same as that of the 
intrinsic curvature radius of the dipolar field lines along 
which they are constrained to move. However, our model 
shows the assumption of identifying the curvature radius 
of particle trajectory with the curvature radius of field 
lines can not be extended to the radio emission region. 

Though it is difficult to find a simple expression for the 
curvature radius to illustrate the effects of rotation based 
on Eq. p3p . we can find an approximate expression 

v 2 
\a\ 

for a particle moving on a field line very close to the mag- 
netic axis. The parameter p m designates the approximate 
expression for the curvature radius when the rotation is 
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Fig. 6. The curvature radius of particle trajectory at different r e of the field lines which lie in the meridional plane (<f) = 0°). 
In each panel, the dotted curve represent the stationary case while the continuous curve represent the rotating. Chosen a = 90° 
and P = 1 sec. 
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Fig. 7. The curvature radius of particle trajectory at different r e of the field lines which lie in the leading (<f> = —90°) and 
trailing sides (cj> = 90°). In each panel, the continuous curve represent the stationary case {cj> = ±90°), while in the rotating case 
the dashed line curve represent the trailing (0 = —90°) side and the dotted one the leading side (</> = 90°) . Chosen a = 90° 
and P = 1 sec. 



taken into account. Consider the approximate forms of 
Eqs. (EM an d CO}: 

d9' d9 cos 6 r 

—— f» COS 0— PS — ■ , 

at dt 2 ^Jr r e 

and 

sin 9' ?s sin a . 

The term dO' /dt k as r c > r for the field lines which 
are close to the magnetic axis. Next by taking v = c(3, 
a = dv/dt and d(f) p /dt as f2 at <fi « 0°, we obtain 

Pm w [4 + (csc 2 a - 1) sin 2 (a„o t)] ~' /2 ■ (15) 

For a = 90° we find p m as tl/2, which is in agreement 
with our earlier estimate in the 2D case (TG05). 



A general expression relating the radius of curvature 
of trajectory of particles moving on a rotating field line: 



^=«( 1+ t)/(v i+ f) ^ 

where pdipoic is the radius of curvature of non-rotating 
dipolc filed line, and the details are given in Appendix-E. 
The expression for Pdipoic is given by Eq. (4) in G04. For 
the altitude very close to the stellar surface both the terms 
v rot and a ro t become negligible in comparison to vq and do- 
But for higher altitudes w ro t and a rot become significant, 
and hence p ro t differs from pdipoic as indicated by Figs. [6] 

km 

Gupta and Gangadhara (2003) by analyzing the pro- 
files of a set of six pulsars indicated that the conal emis- 
sions are coming from the field lines whose foot points 
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Fig. 8. The characteristic frequency (v) and power (P) emitted by a charged particle moving along the field lines with different 
r e , in both rotating and non-rotating cases, are plotted as functions of r/n- In each panel, the rotating case is plotted with 
a solid line while the static case with a dashed. The power emitted is normalized with P max = 6.5 x 10 -26 ergs/sec. Chosen 
a = 90°, 7 = 350, P = 1 sec and (f> = 0°. 



on the polar cap are within the range of approximately 
0.2 to 0.7 Si,. The polar cap radius Si = v/rf /Vl, where 
r s ?s 10 Km is the neutron star radius. For typical classical 
pulsars the value of Si lies in the range of 100 to 200 m, 
and the conal rings may be separated by 10 to 20 m on the 
polar cap (Gupta & Gangadhara 2003). So, in the picture 
of nested conal structure, the central component (core) 
is expected to occupy a region with radius of about 15 to 
20 m around the magnetic axis. We find the constant r e of 
field lines associated with different emission components 
approximately lie in the range of 5 to 70 ri for the conal 
components and greater than 70 ri for the core. 

The characteristic frequency (y) and the power (P) 
emitted by a particle due to curvature radiation (RS75) 
are given by 



47T p 

and 



3 c \p 



(17) 



(18) 



In Fig. [5J we plotted v and P in both the rotating (continu- 
ous lines) and the non-rotating (dashed lines) cases. When 



the rotation is considered we find significantly enhanced 
emission compared to the non-rotating case. In the case 
of emissions from the field lines close to the magnetic axis 
the disparity in the values of frequency and power emitted 
between the rotating and non-rotating cases is huge. This 
again emphasizes the importance of including the rotation 
on pulsar radio emission. 

Many normal pulsars have the average pulse profiles 
with very prominent core component. If the core emission 
is believed to occur from the field lines, which are close to 
the magnetic axis then it becomes difficult to account for 
the strong core due to the intrinsic curvature emission, 
which do not invoke rotation. This is because the field 
lines arc almost straight (pi ~ oo), and hence the intrinsic 
curvature emission is expected to become inefficient. 

Other aspect of p m to be noted is its dependence on 
the inclination angle a and the pulsar period P, and it 
can be easily understood from Eq. (fT5)) . We find that p m 
decreases with the increase of a (see, Fig. [9j: minimum for 
orthogonal rotators but maximum in nearly aligned ones. 
But p m increases with the increase of pulsar period P (see, 
Fig. ITU)) : larger for classical pulsars than in millisecond 
pulsars. However, the pi does not have any such variations. 
In literature, p m of the particle trajectory has been often 



8 



R. M. C. Thomas and R. T. Gangadhara: Radio Emission in Pulsars 




0.02 0.04 0.06 0.08 0.1 



r/r L 

Fig. 9. The curvature radius vs r/rL at different angles of 
inclination a for particles following the field line with r c = 
15rL- Chosen = 0° and P = 1 sec. 




r / 10 cm 

Fig. 10. The curvature radius vs r/r^ at different pulsar ro- 
tation periods P based on Eq. (|13p . Chosen a = 90°, = 0° 
and r e = 15. 



taken to be of the order of 10 9 cm (e.g., Zhang & Cheng 
1995; Harko & Cheng 2002; Gil et al. 2004) or assumed 
to be same as pi. We emphasis that by making such an 
approximation we miss out an important physics behind 
the influence of rotation. 

Asymmetry in the curvature radius of particle trajec- 
tories on leading and trailing sides with respect to the 
magnetic axis is quite evident from Fig. 5. The leading 
side trajectories are found to have smaller curvature radius 
than those on the trailing side. These differences arc fur- 
ther reflected in the generation of asymmetric profiles in 
such way that the leading side components become more 
intense compared to those on the trailing side in a pulse 
profile. This feature was first remarked in BCW91 that 
out of the 88 cone-dominated pulsars 53 pulsars have the 
leading component stronger while only 24 pulsars have 
the trailing component stronger. One can intuitively un- 
derstand this as a consequence of the fact that rotation 
increases the curvature of particle trajectory on the lead- 
ing side but reduces on the trailing side. Thus owing to 
rotation, particle trajectory becomes more curved on the 



leading side but stretched out on the trailing side, as in- 
dicated by Fig. 2. From Eq. (fl4|) we observe though \v\ 
of particle remain much closer to the speed of light c, the 
acceleration \a\ will be much less on trailing side than 
on leading as indicated by Fig. [7] Recently, Pcyman and 
Gangadhara (2002) have analyzed the profiles of 24 pul- 
sars, and reported that in 19 pulsars the leading side com- 
ponents are broader than their trailing side counter parts. 
It is in clear agreement with Fig. 4 that due to rotation 
the relativistically beamed emissions on the leading side 
are spread out in angle compared to the trailing side. 

The characteristic frequency of the curvature emission 
is inversely proportional to the curvature radius of par- 
ticle trajectory (RS75). Hence from Figs. [6] and [7] we ex- 
pect that the high frequency radiation arises from lower 
altitudes than the lower frequency one. This implies the 
existence of radius-to-frcquency (RFM) mapping in the 
case of emission from field lines close to magnetic axis. 
If core emission occurs from these field lines then we ex- 
pect that the core emission heights also follow the RFM. 
Since the curvature radius increases steeply at lower alti- 
tudes (r/r^ < 0.01), most of the high frequency radiation 
is expected to be generated over a narrow range of alti- 
tude. But at higher altitudes curvature radius increases 
less steadily with respect to r. Hence the lower frequency 
emission region is found to be distributed over much a 
larger range of altitude than the high frequency radia- 
tion. It is in clear agreement with the conclusion drawn 
by Hocnsbroech and Xilourics (1997) that the RFM is ev- 
ident at cm-wavelengths but saturates progressively with 
the increase of frequency. 

Since the millisecond pulsar rotates much faster than 
the normal ones it is natural to expect that the effects 
of rotation to be more prominent in them. An easily de- 
tectable and a direct prominent effect of rotation in mil- 
lisecond pulsars is the larger aberration-retardation effect 
in comparison to normal pulsars. For example, one can 
look at our recent work on a nearby and strong millisec- 
ond pulsar PSR J0437-4715 (Gangadhara and Thomas), 
where we find that the aberration-retardation phase shift 
is as high as 23°, a value much larger than the one ever 
detected in normal pulsars. Further example is the work 
of Johnston & Weisberg (2006) on PSR J1015-5719, which 
has a spin period of 139.9 msec. They have detected the 
aberration-retardation phase shifts, which are much larger 
than the normal pulsars (Gangadhara & Gupta 2001; 
Gupta & Gangadhara 2003). 

3.3. Numerical solution 

To estimate the accuracy of our analytical solution r = 
ro + ri , wc numerically solved Eq. ([6]) using Eq. Q . The 
Eq. was split into two first order equations given by 

p = f 

and p = f (r, r) , (19) 
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Fig. 11. The analytical solution (r) up to the first order in e of the radial equation of motion is compared with the corresponding 
precise numerical solution (r„) for different field line constant r e . The error factor <5(%) = (1 — r/r n ) x 100 is plotted with 
respect to the fractional radial distance r/n.. Chosen a — 90° and pulsar period P = 1 sec. 



where p is a dummy variable and 



f(r, f) = - 



20 



2 rf 2 



{dnlJdt)rr 2 



c 2 (l-n 2 m r 2 /c 2 ) 



+ n 2 r 



(20) 



By choosing a = 90° and = 0°, we numerically solved 
the above system of coupled first order equations using 
Matlab ODE-solvers. They have been set as a bench mark 
for testing the accuracy of the perturbative solution. To 
compare the numerical solution (r n ) with the analytical 
solution (r), we have plotted the error factor 5(%) = (1 — 
r/r n ) x 100 with respect to the fractional radial distance 
(r/r^) in Fig. [TT] It shows error factor 6(%) is less than 
0.025% in the radio emission region. 

4. Conclusion 

By taking into account of pulsar rotation, we have devel- 
oped a 3D model for the motion of charged particles in the 
radio emission region of pulsar magnetosphere. We have 
solved the equation of motion by perturbative method and 
obtained the trajectory of particles. For an observer in the 
laboratory frame, trajectories are found to have clear de- 
viations from the geometry of dipolar field lines. Due to 
rotation, the trajectories become more curved on the lead- 
ing side of magnetic axis than on the trailing side. Particles 
emit curvature radiation due to the curvature introduced 
into their trajectory by both the intrinsic curvature of 
dipolar field lines and the rotation. The emissions on the 
leading side are expected to dominate over that on the 
trailing side. We expect that the emissions close to the 
magnetic axis (core) includes a major contribution from 
the rotation induced curvature emission. Our model pre- 
dicts radius-to-frequency mapping for both the core and 
the conal emissions. 
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Appendix A: The small angle approximations 

For analytical simplicity we introduce the following small angle 
relations: 

9' — 9d + a and 
cos 6' — cos(#d + a) 

— cos a cos 9d — sin a sin 9d, (A.l) 

where 9d is the projected angle corresponds to 9 on the merid- 
ional plane. Using the above equations and Eq. ()C-15[1 and 
taking that both 9 and 9d are small we find 



sin Ud ~ cos (p sin ( 
or 9d ~ 9 cos cj> . 



The colatitude with respect to magnetic axis 

9 = sin -1 (i/r/r e ) 
or 9 ps \J r/r e . 

The rate of change of colatitude 
d9 1 1 dr 



(A.2) 
(A.3) 



(A.4) 
(A.5) 



dt 2 r (j- e _ r ) 
1 1 dr 



dt 



2 y/f r c dt 

We also rewrite sin 2 9' as 

sin 2 0' = sin 2 (a + 9d) 
1 



(A.6) 



1 — cos(2q) cos(29d 



+ sin(2«) sin(26>d) 



2 cqs(2q) 

1 +°d . 2 



f- cos 2 ^ (1 - cot 2 a) . (A.7) 



Appendix B: The expression for Magnetic 
Lorentz force 

Since the field line tangent curves off from the direction of 
particle velocity (v\\), an offset angle introduced between v\\ 
and the field line tangent b. This produces a non-vanishing 
magnetic Lorentz force Fb, which in tern constrains particle 
to the field line. The expression for Fb in a 2D case has been 
derived by Gangadhara (1996). Here we are extending to find 
the relevant expression for Fb in the 3D case. 

The rotation velocity of a particle moving along a field line 
is given by 



v TO t{r) = CI x r = rfisin#'e 



(B-l) 



where 9' is the angle between rotation axis Cl and the radial 
vector r, and e is a unit vector in the direction of rotation 
velocity. If O is the angle between the field line tangent 6 and 
e, then 

e = cos Gen + sin Qe± 



where the unit vectors e|| and e± are parallel and perpendicular 
to b, respectively. The expression for has been derived in 
G05. 

The component of rotation velocity of particle perpendic- 
ular to the field line at a distance r is given by 



v p = Vrot(r) sinGex, 



(B-2) 



while the rotation velocity of field line at r + 8r is given by 



Vf = v rot (r + Sr) sin Qe± 



(B-3) 



in the laboratory frame. Since the magnetic field is corotating 
with the pulsar different points on the field lines have different 
velocities for an observer in the laboratory frame. The particle 
at radial distance r sees that the field line bends away from 
its velocity direction. Hence the relative velocity between the 
particle and the field line is given by 



Vrel 



Vf = SrQ sin 9 sin O ej 



(B-4) 



where we have taken that both sin O and sin 9' remain roughly 
constant over an incremental distance Sr. 

Since the gyration frequencies are very large compared to 
the radio wave frequencies the equation of motion Eq. (JTJ can 
be averaged over period of gyration particle (BCW91). Hence 
the effective magnetic Lorentz force is given by 



F B = - Orel x B) 
c 



(B-5) 



Appendix C: 
motion 



The Derivation of the Equation of 



In the superstrong magnetic field of pulsar, charged particles 
lose their initial velocity component perpendicular to the mag- 
netic field via synchrotron radiation. Hence the particle motion 
becomes effectively one dimensional along the magnetic field 
lines. We derive and solve the equation of motion of particles 
moving along the dipolar field lines. 



Cl. Particle motion along rotating magnetic field 

Consider a particle or a bunch moving along the magnetic field 
lines of an inclined and rotating dipole. For an observer in the 
co-rotating frame, the velocity of particle is given by 



/3\\cbt 



(Cl) 



where /3ii is a free parameter, reflecting the fraction of particle 
velocity along field line, b t = bt/\bt\, and b t = dr/ 89 is the tan- 
gent to magnetic field line of an inclined dipole (Gangadhara 
2004; hereafter G04). The position vector of the particle is 
given by r = re T , where r = r e sin 2 9, r e is the field line con- 
stant, 9 the magnetic colatitude, and e r the unit vector (see 
Eq.lC5l below). 

Next for an observer, who is at rest in the laboratory frame, 
the velocity of particle is given by 



«lab 



vu + S~l X r 



where f2 is the pulsar angular velocity. 
By differentiating 



(C-2) 



(C-3) 
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with respect to time in the laboratory frame, we obtain 
dr „ d 



. „,dd> p „ 

viab = -rr e r + r —— e g + r sm 8 —-— e 

at at at 



By taking the dot-product of with Eqs. ()C-2|I and ()C-4[) . 
and equating the right hand sides, we obtain 



(C-4) 



where (e r , eg, e$) axe the basis vectors of spherical polar coor- 
dinate system centered on the rotation axis Cl. Since the equa- 
tion of motion, i.e, Eq. (JXJ) can be averaged over a gyration 
period of particle, the radial distance r is effectively designat- 
ing the guiding center location. 

The magnetic colatitude 8 is measured in the coordinate 
frame with z axis taken along magnetic axis rh while the primed 
variable (colatitude) 8' is measured in the coordinate frame 
with z axis taken along Q. The basis vectors and their deriva- 
tives with respect to time are given by 



e$ 

de r 
~dt 
dee 
~~dt 

de^ 
Ttf 



sin 8' (cos <f>p x + sin 4>p y) + 
cos 8' z , 

cos 0' (cos <fip x + sin <f> p y) — 
sin 6' z , 

— sin cf> p x + cos 4> p y , 

dO' A d(j>p ,„ 

—rr eg H rr sm # e <t> > 

dt at 



cos w , 



- — Me f + 



cos 8' eg) 



C.2. Particle phase 

The particle phase <j) p with respect to 
system-XYZ (see Fig. 1) is given by 

4>p = Qt ± A(f> , 



(C-5) 

(C-6) 
(C-7) 

(C-8) 
(C-9) 
(C-10) 

fixed co-ordinate 
(OH) 



where is the phase of particle measured with respect to the 
phase of magnetic axis rh on the equatorial plane. The positive 
sign is for the leading side and negative for the trailing. The 
vector rh can be defined as 

rh = sin a cos(fit) x + 

sin a sin(ftt) y + cos a z . (C-12) 

If 8 is the magnetic colatitude of particle then we have 

e r • rh = cos 8 
— cos a cos 8 + sin a sin 8' cos(A0) . 

(CMS) 

By taking dot product of the vectors Q and e r , we obtain 
f2 • e r = cos0' , (C-14) 
and the expression for cos#' is given recently in G05: 
cos 9' = cos 8 cos a — sin 8 sin a cos . (C- 15) 

Using Eqs. (|C-13|I and (|C-15|I . we find an expression for A<f> : 
A<f> = 

/ cos 8 sin a + cos a cos sin t 



(CM6) 



rsm^ = (nxr).e, + -6 t 

Next, by simplify it, we obtain 
cos O 



ft +/3||C- 



(C-17) 



(C-18) 



where f2 p = dcj> P /dt and cosG = 6t • e^>. The expression for 
cos 6 is given in G05 (see eq. 5). 

C.3. Finding the equation of motion 

The momentum of particle is given by p lab = mv\^, where 
m = mo 7, and the quantities m and mo are the relativistic 
and the rest mass of particle, respectively. The Lorentz factor 
7 is given by 



7=1-^ 



-1/2 



(CM9) 



By substituting for viah from Eq. (|C-4[l . we obtain 

r dO' 
c dt 



-1/2 



(C-20) 



where f is the time derivative of r with respect to time. Next, 
by substituting for p lab into Eq. |TJ and differentiating with 
respect to time, we obtain 



d 
Tit 



(4) 



mr Sip sin 2 8' 



+ 



d8^_ 



dO' 
It 



dr<W_ 
''~dl~dt 



mrQp sin 8' cos 8' 



m f2. 



d 



r sin 8 ) 



+ 



d ( 



dt V 
= qE + F 



dt 

I m r Q p sin 8' 



(C-21) 



For the convenience of understanding the motion, we resolve 
the equation into radial and non-radial parts, i.e., along e r 
and in the plane: (eg, e^>). Hence, by separating the terms of 
Eq. ()C-21|) into radial and non-radial parts, we obtain 



d 

Tit 



dr\ 

m Tt) 



dO' 

Tit 



Far + (qE r + F CI 



and 
F v -- 



FbOcP + qEg 



(C-22) 



(C-23) 



where Fbi and F B s^ = -Fbs + F B <t, axe the radial and 
non-radial components of Fb, respectively. Similarly, E r and 
Eg$ — Eg + E$ are the components of E. The quantities Fbt, 
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.Fee and Fb4> are the components of F b in the directions of 
(e r , eg, e^), respectively. Similarly, E has the components in 
these directions. 

The centrifugal force is given by 



mflpfi x (r x f!) 
mr sin 8' Slpi p , 



(C-24) 



where 



+ cos 



?' eg = cos ( 



, x + sin <pp y. 



Then the component of centrifugal force in the e r direction 
is given by 

Fcr = (Fc ■ £r) e r . 

The quantity f u is the sum of the Coriolis force, the force 
generated due to non-uniform rotation velocity v rot and the 
component of centrifugal force in the eg and directions: 



dtf_ 

dt 



+ m 



drd9^_ 
dt dt 



mrD.p sin 0' cos 8' 



mQ.„ — (rsinf?') + 
dt 



eg + 



d 

Tit 



(C-25) 



Since the detailed study of particle motion is beyond the 
scope of this work, we intend to qualitatively understand the 
constrained motion. The following deductions are based on the 
presumption that the charged particles closely follow the field 
lines over the range of radial distance of our interest, i.e., the 
radio emission region as mentioned earlier. 

Initially, particles can have a velocity component (v±) per- 
pendicular to the magnetic field. In the superstrong magnetic 
field of pulsar, particles lose their perpendicular component of 
motion via synchrotron radiation losses. Therefore, they essen- 
tially have one dimensional motion along the field lines. Since 
particles move along rotating field lines, there will be a rela- 
tive non-zero rotation velocity v Te \ between the particles and 
the field lines. Hence particles try to gyrate due to the action 
of magnetic Lorentz force: Fb = q/c{v Te \ x B), but the field 
lines curves off to the side. The motion of particles in curved 
magnetic field lines is discussed by Jackson (1975). The expres- 
sion for v Te \ is derived in Appendix B. The forces F\j, Fb84> 
and qEg,/, act as constraining forces. They balance in such a 
away that the particle is hardly allowed to deviate away from 
the field line. This bead-on-wire approximation holds till the 
aforementioned force balance is maintained. The bead-on-wire 
approximation is expected to fail near the light cylinder, as the 
inertial forces such as Coriolis force exceeds the Lorentz force. 
Because the rotation speed approaches the speed of light and 
magnetic dipole field becomes weaker near the light cylinder. 

If Ob is the angle between e r and b, then we have 



sin Hb 



(e r x b t ) ■ (e r x b t ) 
^2sin6» 



a/5 + 3cos(2(9) 



(C-26) 



We find sin #13 <C 1 within the region of 10% tl in the case of 
normal pulsars, but 30% 7"l in millisecond pulsars, therefore, 
we infer that \F B r\ < l-FfleJ. 



By substituting for F CI into Eq. ()C-22[) . we obtain the equa- 
tion of motion in radial direction: 
d ( dr s 



d ( dr 
-r ni- 
di \ dt 



~dt 

= qE T + Fi 



i 2 9' ttl 



(C-27) 



The component E T ~ E\\ is negligible as the component of elec- 
tric field parallel to the field lines is screened by the (e + , e~) 
plasma. Since 8b ~ 0°, the force term Fbt becomes negligibly 
small compared to the other terms in Eq. (|C-2Tfl . Hence the 
equation of motion in the radial direction reduces to 



d ( dr\ 

— m — — m r 

dt V dt) 



dff_ 
~dt 



(C-28) 

In the limiting case of 9 — 0° and d9/dt — 0, i.e., for a 
particle following the magnetic axis of an orthogonal rotator 
(a = 90°), we can show that the Eq. ()C-28|) reduces to the 2D 
equation of motion derived by Gangadhara (1996). 

Appendix D: The Solution of the 
Equation C^Igji 

Here we develop a perturbative method for solving Eq. (|6]). We 
consider e = r-^jr^ as perturbative parameter. 

The series expansions of r and f in terms of e are 



r = ro + e n + e r2 + • 
and 

f = ro + e ri + e r-2 + 



(D-l) 



(D-2) 



Thus by substituting Eqs. ()D-1|) and (|D-2f) into Eq. 0, and 
series expanding it we obtain 

tirn = ^mO + 6 ^ml + ^ ^m2 + ' ' ' (D-3) 

and 

d£L _ 
dt 

"mO 



(D-4) 



where the perturbation expansion coefficients, which follow 
from Eqs. ©, (|A~Tj) and (fA~6)) . are 



o 2 


= 


• 2 

sin a , 








n 2 

"ml 


= Ri 


' l +R 2 




ro r L 


n 2 


= 




n 2 , 

ml 


= Ri 


( 2r r 




\ r r L 


and 






Ri ~- 




- cos(2(j})] 



r\ r L 



+ R^ 

rL 



(D-5) 



R2 = -r 



1 + 3 cos(2a) 

2 cos (2 <f)) sin 2 



(D-6) 
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After substituting Eqs. (|D-1|I - (|D-6|) into Eq. (®, the coef- 
ficients of e give the zeroth order equation: 



ro + 



and those of e give the first order equation: 
4 fi 2 l0 r r 



(D-7) 



ri + 



+ 



+ 



^(l-fi^r 2 /^) 
4 fit 



n 



l m0 r r 



ri 



c 4 (l-fiL,r 7c 2 ) 2 



2 2r2 ml r r r2 2 nl ror 



^(l-^o'-oVc 2 ) 
"c 4 (l-^o 2 /c 2 ) 2 ■ 



(D- 



D.l. The solution 



We employed the method, which has been used in our 2D model 
(TG05), to find the analytical solution to zeroth order equa- 
tion. Thus the solution of the zeroth order equation is given 
by 



r 



-cn(A - fi m0 t) , 



(D-9) 



where cn(z) is the Jacobian Elliptical cosine function 
(Abramowicz & Stegun 1972) and 



'ill 



A = 



dC 



,Jl - Psin 2 C 

The functions k and (fro are given by 
1 



k z = 

1 - sf, 

and 

4>o — arccos 



so 



(1 - sl)W 



(D-10) 



(D-ll) 



where sq — s| ro=r , , so = ds/dt\t=o, and s — £l m oro/c. The 
parameter n ~ 10 6 cm is the initial radial location of particle 
at time t = 0. The expression given by Eq. (|D-9[l is an exact 
solution for particles moving on the magnetic axis, and it is 
applicable to both the normal and the millisecond pulsars. We 
substitute the approximate form of Eq. (|D-9[) into Eq. (|D-8[I 
to find a solution to the first order equation. Since k ~ 0, we 
can approximate ro and obtain 



ro 



-sin(fi m ot) 



(D-12) 



We seek for a solution for the first order equation which is valid 
over a distance 10% rn from the neutron star surface. Since 
Q m t <C 1 over this range of altitude, we can take sin(fi m ot) ss 
ttmot and 



r 2 O 2 
'0 S ! md 



Using the approximations (|D-12f) and ()D-13[) we can rewrite 
Eq. (EDI : 



(D-14) 



r\ + AQtth +Qi(l + 4Qit 2 )n 



= P 2 t 2 + P 4 t* 



where 



1 — f^mO • 

2c 2 (i?ifi m0 2 



P 2 = - 
Pi = - 



R2) 



2C 2 -R2^m() 2 



The homogeneous part of Eq. (|D-14|) is given by 
/i + 4Q 1 t/ l + gi(l+4Qit 2 ) h = 0, 



(D-15) 



(D-16) 



where h is the contribution to n from homogeneous part. By 
defining 



ft = /io(*) exp[— Qi t 2 ] 
we can reduce Eq. (|D-16|) : 
ho — Qi ho . 

The solution of Eq. (|D-18|) is given by 
ho = Ci sinh(-\/(27f) + C2 cosh( — yQit) 



(D-17) 
(D-18) 

(D-19) 



where Ci and C2 are arbitrary constants. We take the initial 
condition that ho = and ho = « p o and -u^o <S c at time t = 0. 
Thus we have G\ = « P o/\/Qi an d C 2 = 0, and the solution /io 
can be written 



2 



r ( 



exp 



exp 



-y/Qit 



where Since exp [\/Qi i\ and exp [— \JQ\ t\ independently sat- 
isfy Eq. (|D-18[) . we can write the solution for ri as 



ri = 
where 

yi = exp 
j/2 = - exp 



2/2 k 



dt + 2/2 



2/1 « 



dt 



(D-20) 



?ii 2 + 



(D-21) 



d 2/i d 2/2 
to = 2/2-7- - 2/1- 



d£ dt 

is the Wronskian (e.g., Zwillinger 1989), and k = P 2 t 2 + P 4 t 4 . 
We numerically evaluated the integrals in Eq. (|D-20|) using 
Mathematica. 



Appendix-E: The relation between pdipoic an d Prot 
By finding the derivative of i>i a b, we can write acceleration as: 

d Vl a b 
a = ; . 

dt 

Hence, following Eq. (|13[) . we find 

1 2 



(E-l) 



(D-13) Prot = 



Vlab 



(E-2) 
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where 

|v lab | 2 =vl+v* ot (E-3) 
and 

2 dr 2 dff 2 . , 2 

vLt = (r-sin6»') 2 (fi 2 + 20F) (E-4) 

where F = j3u ccos 0/r sin 0'. The expression for a can be ex- 
pressed as 

a = a r e r + ag eg + a# e$. (E-5) 
We can split the components as 

a r = a r o + Hrrj 
019 = ISO + flflr, 

(14, = ai0o + a 0r, (E-6) 
where 



d d0'\ drd6' 

aeo = dt[ r ^t) + dt^t 

7l2 



r .F sin cos 8' 
a^o = F— (rsinfl') + 

-|(r-Fsm0') (E-7) 

and 

a„ = -r sinV(ft 2 +2QF) 
a gr = -r sin6»'cos6>'(fi 2 +2QF) 

We write the term a = */af + + a? = ya§~~i~ a rot an d 



2f7^(rsin6''). (E- 
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a = a r0 + a go + a^ 

a 2 ot = 2(a r oa rr + a^oi^r + agn.ag r ) 

+ a 2 r + a 2 .,. + a\ r . (E-9) 

In the absence of rotation, the radius of curvature of parti- 
cle trajectory becomes same as the radius of curvature of field 
line along which it is constrained to move. We can also find 
that the terms v 2 ot and a 2 ot vanish in the absence of rotation 
because of the Q term. Thus we can write 

Pdipole = — • (E-10) 

After a few algebraic manipulations with the above said terms 
and using Eq. (|E-2|I we can write the relation between p rot and 

Pdipole &S 

/ l + v 2 ot /v 2 \ 

Prot = Pdipole . == f ■ (E-ll) 



V 1 + a rot/ao 



